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ABSTRACT 

\ 

The  formulation  of  optical  signature  features  useful  in 
discrimination  studies  and  decoy  design,  and  a  method  for  evalu¬ 
ating  a  decoy's  effectiveness  are  presented  in  this  report.  These 
features  are  based  on  an  optical  signature  model  which  incorporates 
a  second-order  trend  with  a  sinusoidal-type  variation  normally 
observed  in  optical  signatures.  The  decoy  evaluation  procedure 
presented  here  provides  the  decoy  designer  with  an  analysis  of 
those  features  and  wavebands  of  the  decoy  signature  which  may 
need  modification  to  achieve  a  prescribed  match  with  an  RV. 
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The  purpose  of  this  report  is  to  introduce  and  document 
results  obtained  in  the  development  and  use  of  features  derived 
from  optical  signatures.  The  primary  motivation  for  and  appli¬ 
cation  of  these  features  has  been  in  the  development  of  discri¬ 
mination  algorithms  and  in  the  requirements  for  decoys.  The 
latter  application  has  resulted  in  the  development  of  a  computer 
program  which  provides  an  optical  decoy  designer  with  an  evalua¬ 
tion  of  a  candidate  decoy  by  means  of  a  detailed  analysis  of 
the  first-  and  second-moment  statistics  of  its  features. 

The  formulation  of  these  features  and  the  details  of  the 
decoy  evaluation  procedure  are  the  subject  of  this  report;  a 
subsequent  report  will  provide  documentation  of  the  computer 
program*.  The  use  of  these  features  in  discrimination  studies, 
as  well  as  in  decoy  design,  together  with  related  topics  will 
be  presented  when  they  become  available. 


*  B.  J.  Burdick  and  M.  R.  Colaneri,  "BODE;  Bhattacharyya 
Optical  Decoy  Evaluation,"  M.I.T.  Lincoln  Laboratory 
(to  hr  nublished) . 
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II.  Signature  Model  and  Features 

As  an  aid  to  decoy  design,  we  have  developed  an  optical 
signature  model  which  incorporates  the  major  features  observable 
by  a  multiband  optical  sensor.  These  signature  features  are 
easily  calculated  and  are  traceable  to  the  thermal  and  dynamic 
properties  of  the  targets  so  that  any  decoy  modifications 
required  to  achieve  a  prescribed  match  with  RV  signature 
statistics  can  be  easily  ascertained. 

The  signature  model  for  a  given  waveband,  AX,  is  written 
in  terms  of  the  features  as 

“aX^"^  “  ^Ax[^  ”®AX  "^Sx  ^  “aX^"^] 

n=l, 2 , . . . ,N 

li=(n-nQ)/nQ  ,  nQ=(N+l)/2 

where  n  is  the  observation  number,  N  is  the  total  niimber  of  obser¬ 
vations,  n  is  a  normalized  observation  number,  and  nQ  is  the  midpoint 
of  the  observation  interval.  The  first  three  features  appear  ex¬ 
plicitly  in  the  above  equation: 

=  temporal  mean!*  of  (n) 

=  linear  trend  of 

=  quadratic  trend  of 

The  remaining  two  features  are  implicitly  contained  in  H. . (n) , 

A 

which  corresponds  to  the  signature  after  trend  removal,  and  are 
*  To  be  precise,  A^j^  is  the  temporal  mean  only  when 
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2 

cf.,  =  variance  of  H.,(n) 

AA  AA 

^(f)  =  temporal  autocovariance  of  H^^(n) 

A=l,2, . . . ,N-1 

The  assumption  has  been  made  that  the  temporal  autocovariance 
function  need  be  calculated  for  one  band  only  (note  that  the 
subscript  AA  is  absent  from  pCJl)).  This  is  justified  for  two 
reasons:  (1)  it  has  been  observed  that,  within  the  allowable 
statistical  sampling  errors,  ^'(f)  is  the  same  for  all  bands, 
and  (2)  it  reduces  the  number  of  features  and  thus  simplifies 
the  task  of  decoy  design. 

These  five  features,  then,  are  used  to  describe  the 
major  characteristics  of  optical  signatures:  the  features 
^AA'  ®AA'  ^'AA  slowly  varying  trend  and  the 

features  and  "^(Z)  describe  the  remaining,  faster  varying, 

sinusoidal- type  variation  normally  observed  in  optical  signatures. 

The  trend  features  are  obtained  by  means  of  a  least- 
squares  fit  to  the  signature.  That  is,  the  residual 

N  2 

^  [“aa  “^aa"^aa"^^'^aa] 

is  minimized  with  respect  to  a^^^  =  ®AA' 

°AA  “  ^AA  ^AA*  Defining 
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(2) 


n=l 
N 

a*  H^^(n)  ,  (3) 

n=l 

and  noting  that 

<N°>  =  N 

<N^>  =  0 

<N^>  =  N(N-1)/3(N+1) 

<N^>  =  0 

<N^>  =  N(N-l) (3N^-7)/15(N+l) ^ 

We  find  that  the  minimization  requirements  yield  the  following 

=*.  =  <H^“> 

=#>  +  b^^<N2>  +  =  <H^^> 

^  ^  *  ^®AX^  • 

The  solution  of  these  equations  is 


equations 


3R 


da 


AX 


aR 


AX 

dR 


AX 


4 


‘AX 


^AX  = 


'AX 


<Ha?><N^>-<H.2><n2> 

<A<nS-<n^^ 


<"w> 

<n2> 

<ha;><n">  -<h,;><kQ> 

<N°>  <N^>  -  <N^>^ 


so  that 


‘AX 


=  a 


B 


AX 


'AX 


AX 

°4X^“4X 


The  features  related  to  the  sinusoidal-type  variation 
are  obtained  from  the  signature  after  trend  removal, 


“aX^"^  “  {^AX^"^  ■  ^Ax[^  "®AX  "  ^Ax]}/^ 
and  are  given  by 

^AX  ^  N  £  “aX^"^ 


AX 


N 


n=l 

N-H 


cr  2 
AX 


(6) 


(7) 


”aX^"^  H^j^(n+a)/a 

n=l 

£*1,2,. . . ,N-1 


For  the  case  when  H^j^(n)  is  Gaussian,  these  features  provide 
a  complete  description. 


* 


To  gain  a  better  understanding  of  these  features  consider 


the  following  excunple: 

(n)  ~ sin(na)At+9)  (example) 


where  At  is  the  time  spacing  between  observations,  f  is  a  random 
phase  having  a  probability  density 


p(<P) 


{l/2rr,  0<9<2n 
\0,  otherwise 


and  0)  is  a  random  frequency  having  a  probability  density 

p(«o)  .  (For  a  spin-stabilized  target  which  is  axially  symmetric, 

(1)  would  correspond  to  the  precession  frequency  and  would 

be  related  to  the  coning  angle.)  The  features  for  this  signature 

are  then  given  by 

~  2  ^  2 

°AA  *  sin^CnojAt+V*) 

n*l 


and 


2  )  .  sin  (a)At-i-2y)  -  sin  [(2N-H)  coAt+2y] 
AX  I  2Nsina)At 


Pd) 


N-£ 


E 


sin  (nctiht+V) 


sin  [(n+i)  a)At+ 


o  ^ 

I  cosicoAt  +  sin  [(i-fl)toAt4-2<i>]  -  sinf(2N-JL+l)(oAt+2Vj 
2  I  N  2Nsina)At 

AX 


(Note  that  in  the  limit  as  N-»oo ;  and  p  (i) -*  cosiwAt. ) 

The  ensemble  averages  of  these  features  are 
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[ 


=  a 


AX 


and 


Efe(il)|  =  /'“  jC^!  P(^)  P«^)  p(a))dlPda) 

I  i  •'(0=-,  00  •'<*’= 0 

=  — 5^  /  p(a))  cosilcoAt  do) 

N  J~ao 

N-X,  „  r  /*°°  /  V  i)la)At  ,  I 

=  -jj— xRe  J  p(a))  e  dcoj 


That  is,  the  ensemble  average  of  "^(S-)  is  the  cosine  transform 
(or  the  real  part  of  the  characteristic  function)  of  the  fre¬ 
quency  distribution  p(a))  .  This  means  that  E|p(il)|^  contains  all 
the  information  concerning  the  frequency  distribution  of  il(n)  . 

For  example,  if  p(a))  is  Gaussian-completely  characterized  by 

2 

a  mean,  m^,  and  a  variance,  a^: 


p(a)) 


(“-"‘a,) 


2 


/o 


2 

CO 


then 


N-jl  _  r  iilm  At-(Jla  At)  ^/2l 
-jj—  X  Re  I  e  CO  CO  '  '  J 

N-S,  -(Ha  At)  ^/2  .  .  . 

—rr-  e  CO  '  cosHm  At 

N  CO 


which  is  a  damped  cosinusoid  with  frequency  equal  to  the  mean 
of  p(co)  and  a  damping  factor  proportional  to  the  width  of  p(co)  . 
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III.  Decoy  Specification  and  Evaluation 

The  method  used  to  levy  decoy  requirements  and  evaluate 
the  effectiveness  of  a  decoy  against  an  RV  is  based  on  the 
Bhattacharyya  distance?  This  metric  provides  a  measure  of  the 
separability  of  the  RV  and  decoy  in  terms  of  the  ensemble 
means  and  covariances  of  their  features.  It  is  particularly 
useful  in  as  much  as  it  is  directly  related  to  the  probability 
of  error,  is  additive  for  independent  features  and  can  be  ex¬ 
pressed  in  terms  of  differences  primarily  in  the  means  and 
variances  in  each  waveband  and  to  the  correlation  between 
wavebands  for  each  feature.  That  is,  if  (A,B,C,a,p)  represents 
the  Bhattacharyya  distance  for  all  features  in  all  wavebands 
then,  if  it  can  be  assumed  that  the  features  are  independent, 
M(A,B,C,?,?)  =  y(A)  +  y(B)  +  y(C)  +  y  (?)  +  y  (?)  (8) 

where  y(*)  is  the  Bhattacharyya  distance  for  an  individual 
feature.  Further,  as  will  be  shown  later,  each  y(*)  can  be 
separated  into  components  which  are  primarily  related  tc 
differences  in  the  RV  and  decoy  means  and  variances  in  each 
band  and  the  average  correlation  between  bands 

It  is  in  this  manner  that  it  becomes  possible  to  pinpoint  not 
only  which  decoy  features  may  be  contributing  to  a  large 

*  K.  Fukunaga,  Introduction  to  Statistical  Pattern 

Recognition  "(Academic  Press,  New  York  and  London,  1972). 


Bhattacharyya  distance  but  also  whether  this  is  due  to  the  mean, 
variance,  or  band-to-band  correlation  and  further,  in  which  band 
it  is  occurring.  The  Bhattacharyya  formulation  also  makes  it 
possible  to  determine  how  the  decoy  requirements  should  be  levied 
on  the  various  features  and  their  ensemble  means,  variances  and 
correlations  in  the  individual  wavebands.  This  simply  becomes 
a  matter  of  apportioning  the  maximum  total  allowable  Bhattacharyya 
distance  among  the  features  and  then  among  each  of  the  feature's 
first  and  second  moment  statistics. 

The  maximum  allowable  total  Bhattacharyya  distance, 

nicix 

is  related  to  the  total  probability  of  error  achieved  by  the 


decoy,  e,  according  to  the  equation 


(10) 


(which  is  an  exact  equation  when  the  RV  and  decoy  covariance 
matrices  are  identical,  and  a  good  approximation  when  they 
are  not  too  dissimilar) . 

For  example,  if  one  wanted  to  design  a  decoy  to  achieve 
at  least  a  15%  error  rate  with  an  RV,  then  a  total  Bhattacharyya 
distance  of  0.5  or  less  would  be  required  (from  the  equation 
above) .  If  this  were  to  be  apportioned  equally  among  the  five 
features  then  each  feature  would  have  to  have  a  Bhattacharyya 
distance  of  0.1  or  smaller.  If  this  were  further  apportioned 
equally  to  the  first-and  second-order  statistics  of  each  feature 
(i.e.  the  mean  and  the  combined  variance  and  correlation). 
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then  they  would  each  have  to  achieve  a  Bhattacharyya  distance 
of  0.05  or  less.  If  all  the  individual  Bhattacharyya  distances 
were  less  than  these  requirements,  then  the  decoy  would  be 
guaranteed  to  achieve  the  desired  error  rate.  However,  it  may 
still  be  possible  for  a  decoy  to  achieve  the  desired  error 
rate  even  if  some  of  the  Bhattacharyya  distances  are  too 
large-all  that  really  counts  is  the  total  Bhattacharyya  distance 
and  some  components  may  compensate  for  others. 

Let  us  now  show  how  the  Bhattacharyya  distance  for  a 
single  feature,  yC*),  can  be  written  as  a  sum  of  components 
attributable  to  differences  in  the  RV  and  decoy  means  and 
variances  in  each  waveband  and  the  average  correlation  between 
wavebands.  In  general,  the  mean  vector  and  covariance  matrix 
of  a  feature,  X,  measured  in  K  wavebands  are  given  by 
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and  o^j  is  the  covariance  of  the  ith  and  jth  components  of  X: 


a 


2 

ij 


E{ (x^-m^) (Xj-m^) } 


a 


2 

ji 


The  covariance  matrix  can  be  written  in  the  following  form  to 
exhibit  the  explicit  correlation  between  wavebands,  ^ : 


'11 


^12'^11°22 


*^1K®11°KK 


^12‘^ll‘^22  •  •  •  ^Ik'^H^KK 

2 

22 


KK 


where 


p .  .  =  0 .  ./a .  .a  .  . 

^xj  13'  11  33 


(This  p  should  not  be  confused  with  the  feature  ^(S,).)  To 
render  the  remaining  calculations  tractable,  it  is  necessary  to 
make  the  assumption  that  the  correlation  between  band  pairs  is 
the  same  for  all  pairs  so  that  the  p. .s  can  be  replaced  by  their 


average : 


P°l°2 


paiOR 


where 


K  K-1 

P  =  Kfl~l)  ^  Pij 


(11) 


2  2 

We  have  further  simplified  the  notation  by  writing  for 
The  replacement  of  the  individual  band-to-band  correlations  by 
their  average  has  been  found  to  be  a  very  good  approximation 
and  is  justified  on  the  grounds  of  simplifying  the  decoy  design 
procedure. 

The  determinant  and  inverse  of  a  matrix  of  this  form  can 
be  shown  to  be 


12^1  =  •••  a^(l-p)^"^  [l+(K-l)  p] 


(12) 


J+(K-2)p]/aJ  -P/Ol‘^2 


>-l  1 

s 


•  -P/c^i^K 


-P/°l°2 


-P/OiCk  . 


[i+(K-2)p]/02 


[l+(K-2)p]/a^ 


where 


=  (l-p)[l+(K-l)p] 


(13a) 


(13b) 
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The  Bhattacharyya  distance  for  a  given  feature  (under  the 
Gaussian  assumption)  is  a  function  of  the  mean  vectors  and  co- 
variance  matrices  of  the  RV  and  decoy  (DY) 

U(-)  =  J  •  IE®’'!  (14«) 

where  AM  is  the  difference  in  their  mean  vectors; 

am  =  M^^  -  M°^  (14b) 

and  Is  the  average  of  their  covariance  matrices: 


2= 


DY 


(14c) 


With  the  assumption  that  the  correlation  between  band  pairs  is 
the  same  for  all  band  pairs  E  ^sn  be  written  as 


E= 


-2 

°1 


paia2 


(15a) 


-2  . 


pCf^Oj^  •  • 

where  is  the  average  of  the  RV  and  decoy  variances  in  band  k; 

(15b) 


«2  _  i,/„RV.  2  ^  .  /„DY,2 


and  p  is  the  average  band-to-band  correlation  of  the  average  of 
the  RV  and  decoy  covariance  matrices : 

j=i+l  i=l 
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The  determinant  and  inverse  of  are  then  given  by  the 
same  equations  derived  previously  so  that  the  Bhattacharyya 


distance  becomes 


where 


*kli  * 


^COR^*^ 


’‘MNk* 


'-I  ■ '  ‘•K'”?'”!'] 


u  (.)  z=  k  In 


( _ (i-o>'^'^Ch-(k-i)o3 _ I 

(K-H  p"'']  p'*’']) 


where  Amj^  is  the  kth  (waveband)  component  of  AM  and 

=  (1-p)  [l+(K-l)p]  .  (20) 

The  individual  Bhattacharyya  distances  (•)  and  (•) 

k  k 

provide  information  on  mismatches  between  the  decoy  and  the  RV 

due  mainly  to  the  mean  and  variance,  respectively,  in  band  k 

while  contains  information  on  mismatches  attributable  mainly 

to  the  average  band-to-band  correlation.  By  "mainly"  we  mean: 

if  the  only  difference  between  the  decoy  and  the  RV  is  the  mean 

in  band  k,  then  the  only  nonzero  Bhattacharyya  distance  will 

be  (•);  if  the  only  difference  is  the  band-to-band  correla- 

k 

tlon,  then  only  will  be  nonzero;  and  if  the  only  difference 
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T 


is  the  variance  in  band  k,  then  not  only  (•)  but  also  VIqor^*^ 

k 

will  be  nonzero.  That  is,  the  variance  and  correlation  components 
cannot  be  separated  when  there  are  variance  differences  between 
the  decoy  and  the  RV.  An  examination  of  the  equation  for  F  shows 
how  this  occurs.  Consider  the  case  when  the  decoy  band-to-band 
correlation  is  the  saune  as  that  of  one  RV  (p  =p  )  and  the 
decoy  standard  deviations  are  the  seune  as  the  RV  except  for  one 
band  for  i^k).  Then  p  becomes 


where 

r  =  i  [<K-2,  . 


SO  that 


^COR^*^ 


*5 


In 


(21a) 


(21b) 


From  this  equation  we  see  that  will  be  nonzero  when  the 

decoy  standard  deviation  differs  from  the  RV  in  band  k  (even 
though  they  have  the  same  band-to-band  correlation) .  Therefore, 
in  this  case,  the  Bhattacharyya  distance  due  solely  to  a  variance 
difference  in  band  k,  will  be 


VAR,^ 


(•)  =  (•)  + 


VAR 


COR 


(22) 


To  gain  a  clearer  understanding  of  the  dependencies  of  these 
various  Bhattacharyya  distances  on  the  statistics  of  the  RV  and 
decoy  features,  plots  of  their  variations  are  shown  in  Figs.  1-3 
for  the  case  of  three  wavebands. 
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Figure  1  illustrates  the  variation  of  the  Bhattacharyya 


distance  due  to  mean  differences  in  a  single  band,  (•),  with 

""k 

the  mean  differences  normalized  by  the  average  standard  deviation, 
for  various  values  of  the  average  band-to-band  correlation, 
p.  Figure  2  shows  the  variation  of  Bhattacharyya  distances  due 
to  variance  differences  in  a  single  band, 

with  the  ratio  of  standard  deviations,  /a^  or  its  inverse,  for 

various  values  of  the  average  band-to-band  correlation,  p  ,  in 

the  case  of  ( • )  .  Figure  3  shows  the  variation  of  the 

Bhattacharyya  distance  due  to  differences  in  the  average  band- 

to-band  correlation,  Vi^Qp(*)/  with  the  RV  and  decoy  average 

RV  DY 

band-to-band  correlations,  p  and  p  .  This  is  presented  as 
a  function  of  a  parameter  n , 
i-p'”' 

"  “  i-p"''  ' 


so  that  a  single  curve  (approximately)  represents  the  variation 

DY  RV 

of  Vij,Qj^(*)  with  p  for  any  value  of  p  .  Shown  at  the  bottom 

DY 

of  this  figure  are  various  scales  of  p  that  can  be  used  for 

RV 

some  particular  values  of  p 

Given  the  RV  feature  statistics  and  the  desired  error  rate 
to  be  achieved  by  the  decoy,  it  is  possible  to  levy  requirements 
on  the  decoy  feature  statistics  using  these  curves.  Consider 
the  ex£unple  presented  earlier  where  a  total  Bhattacharyya 
distance  of  0.5  was  required  (corresponding  to  a  total  error 


0 


Fig.  1.  Variation  of  the  Bhattacharyya  distance  due  to  a 
difference  in  means  in  band  k  for  several  values  of  aver¬ 
age  band-to-band  correlation. 
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Fig.  2.  Variation  of  Bhattacharyya  distances  due 
to  variance  differences  in  band  k,  with 

shown  for  several  values  of  average  band-to-band 
correlation. 


Bhattacharyya  distance 
irrelation  differences. 


L9 


a 


probability  of  15%).  If  this  is  split  up  equally  eunong  the  five 
features,  then  the  Bhattacharyya  distance  for  each  feature  should 
have  a  value  of  0.1  or  less.  Splitting  this  up  equally  between 
the  first  and  second  moment  statistics  results  in  a  value  of  0.05 
for  the  mean  and  0.05  for  the  variance  plus  correlation  (which 
can  further  be  split  equally  among  the  variance  and  the  correla¬ 
tion  to  give  0.025  for  each).  If  there  are  three  wavebands,  then 
the  mean  and  variance  distances  should  be  apportioned  equally 
among  the  three  bands  so  that  we  arrive  at: 

3c 


If  the  average  correlation  between  bands  is  0.90  for  the  RV  then, 
reading  from  the  plots,  this  means  that  the  decoy  feature  sta¬ 
tistics  must  satisfy: 


mj^  -  0. 14aj^  <  m^^<  in^^+0. 14CTj^ 


0.91of<  of  <1.10of 


DY , 


0.84  0.936 

These  would  be  the  requirements  that  each  of  the  decoy's  features 
in  each  band  would  have  to  meet  in  order  for  the  decoy  to  be 
assured  of  achieving  a  15%  error  rate.  This  is  an  illustrative 
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example  only;  the  precise  decoy  requirements  will,  of  course, 
depend  on  the  RV's  average  band-to-band  correlation  for  each 
feature  and  on  the  total  error  rate  desired  for  the  decoy. 

These  results  apply  only  to  the  features  ®AX'  ^AX' 

o the  Bhattacharyya  distance  and  the  decoy  requirements  for 
the  temporal  autocorrelation  feature, p (H),  are  computed  in  a 
different  manner.  There  are  basically  two  distinct,  and  comple¬ 
mentary,  ways  in  which  the  Bhattacharyya  distance  corresponding 
to  the  RV  and  decoy  'Jf(iJ.)s  can  be  calculated. 

The  first  method  consists  of  computing  the  spectra  of  the 
RV  and  decoy  "pCiDs,  estimating  the  mean  and  standard  deviation 
of  these  spectra : 


"'(0  '  spectrum  mean,  standard  deviation 

DY  DY 

m^  ,  =  decoy  spectrum  mean,  standard  deviation 

and  computing  the  univariate  Bhattacharyya  distance  (which  is 
again  separable  into  mean  and  a  variance  contribution) : 

ti(P)  -  P^«.)  +  (23a) 


=  *5  In 


o2 

0) 


^RV^DY 
a  a 
0)  01 


(23b) 

(23c) 


*  For  the  feature  statistics  should  be  computed  for  ^■naAX 
since  it  is  more  likely  to  be  lognormally  distributed  than  normally 
distributed  (i. e. » negative  values  are  not  allowed). 
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where 


.  RV  DY 

Am  =  -m,, 

0)  0)  CO 


-2  2  ^  i,/„DY.2 


(23d) 

(23e) 


DY  DY 

The  levying  of  decoy  requirements  on  m^^j  and  can  be  achieved 
in  the  same  manner  as  described  for  the  other  features  using  the 
plots  of  Figs.  1  and  2  (with  P=0) .  This  method  is  useful  since 
it  provides  statistics  oh  the  frequency  of  the  sinusoidal-like 
variation  of  the  signature.  However,  it  assumes  that  the  frequency 
distributions  are  unimodal  and  can  thus  be  characterized  by  a 
single  mean  and  standard  deviation. 

The  second  method  is  based  directly  on  the  first-and  second- 
moment  statistics  of  the  sinusoidal-like  signature  remaining 
after  trend  removal,  H^^(n),  and  makes  no  assumptions  about  the 
nature  of  the  frequency  distribution.  If,  after  trend  removal, 
the  signature  is  normalized  to  unit  variance  (by  extracting  the 
feature  then  it  is  completely  characterized  (for  Gaussian 

statistics)  by  its  ensemble  mean  vector  and  covariance  matrix, 
which  are  given  by; 


M  = 


0 


0 
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1 


E{?r(i)i  . 


.  E{^(L)} 


E{p(l)}  1 


E{p(l)  } 


E{p(L) } 


E{^(1)}  1 


(24) 


where  L  is  the  number  of  lags  (which  is  less  than  or  equal  to 
N-1,  where  N  is  the  number  of  observations) .  The  mean  vector  is  zero 
because  the  trend  has  been  removed  and  the  covariance  matrix 
has  the  form  that  it  does  because  of  the  definition  of  'p()l). 

That  is,  an  element  of  21  is  defined  as 

and,  if  we  can  assume  stationarity ,  is  a  function  only  of  the 

difference  i!,=  lm-n|  so  that 


Interchanging  the  sximmation  and  expectation  yields 

=  E{^(]1)}  . 

It  should  be  noted  here  that  we  are  using  a  biased  estimate  for 

9" _ and  9(Jl)  since  the  sums  are  divided  by  N  and  not  by  N-Jl 

(which  is  the  actual  number  of  terms  in  the  sum) .  This  is  done 
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to  ensure  that  the  covariance  matrix,  S  ,  is  positive  definite. 

The  Bhattacharyya  distance  corresponding  to  the  feature 
?(£)  is  then  given  by  (only  the  variance  term  contributes  since 


the  means  are  zero) : 


where 


M(P)  =  %  in  if  I  A/®" 

f=  . 


(25a) 


NOTATION 


waveband 

total  number  of  wavebands 

waveband 

observation 

frequency 


decoy 

reentry  vehicle 


normalized  or  based  on  normalized  signature 
average  (over  wavebands  or  of  RV  and  DY) 


lag 

total  number  of  lags 
observation 

total  number  of  observations 
midpoint  of  observation  interval 
normalized  observation  number 


NOTATION  (cont.) 


Operations,  Functions 


<  > 

suimnation  over  observations  (temporal  sum) 

E{*} 

ensemble  average 

transpose 

1*1 

determinant 

inverse 

p  ( • ) 

probability  density 

Variables 

temporal  mean  feature  (of  H^j^(n)) 

^AX 

=^X 

®AX 

linear  trend  feature  (of  H^^(n)) 

^AX 

"^^AX^AX 

^AX 

quadratic  trend  feature  (of  (n) ) 

°AX 

“^AX^AX 

e 

probability  of  error 

"sx'"’ 

optical  signature 

»4X<"> 

normalized  optical  signature  (zero  mean, 
unit  variance) 

temporal  moment  of  H^^(n) 

M 

ensemble  mean  vector  of  a  feature 
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NOTATION  (cont.) 


Variables  (cont.) 


am 

difference  in  RV  and  decoy  ensemble  means 

M 

ensemble  mean  of  H..(n) 

Aa 

waveband  component  of  M 

Amj^ 

waveband  component  of  AM 

•  mean  of  spectriam  of  ^(£) 

Am 

0) 

difference  in  RV  and  decoy  spectral  means 

y(-) 

Bhattacharyya  distance 

Bhattacharyya  distance  due  to  mean 
differences 

^VAR^*^ 

Bhattacharyya  distance  due  mainly  to 
variance  differences 

^COR^*^ 

Bhattacharyya  distance  due  mainly  to 
correlation  differences 

^COR  ^  ^ 

Bhattacharyya  distance  due  to  variance 
differences  reflected  in  P^OR^*^ 

^var(-^ 

^VAR^*^  ^COR^*^ 

y 

max 

maximum  allowable  total  Bhattacharyya 
distance 

n 

=  (1-p  )/(l-p  ) 

<N*> 

sum  of  n* 

random  phase 

R 

residual  in  mean-square  minimization 

p(4) 

temporal  autocovariance  feature  (of  (n) 
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NOTATION  (cont.) 


Variables  (cont. ) 


P 

"p 


E 


a 


2 

ij 


-^2 

*^nm 


0) 

X 


X 


k 


correlation  between  wavebands  i  and  j 

average  band-to-band  correlation 

average  of  RV  and  decoy  average  band-to- 
band  correlations 

ensemble  covariance  matrix  of  a  feature 
average  of  RV  and  decoy  covariance  matrices 
ensemble  covariance  matrix  of  H^^(n) 

i , j  element  of  E 

2 

simplified  notation  for  Oj^j^ 

2 

average  of  RV  and  decoy  Oj^  s 

temporal  variance  feature  (of  H^j^{n)) 

m,n^  element  of  2 

=  t  i=lm-n| 

mn  ' 

variance  of  spectrum  of 

average  of  RV  and  decoy  spectral  variances 
standard  deviation  of  exeunple 
time  spacing  between  observations 
frequency 
arbitrary  feature 

waveband  component  of  X 
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